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Abstract 

We present a theory for the spin correlation function of the t-J model in the 
framework of the dynamical mean-field theory. Using this mapping between the 
lattice and a local model we are able to obtain an intuitive expression for the 
non-local spin susceptibility, with the corresponding local correlation function as 
input. The latter is calculated by means of local Goldstone diagrams following 
closely the procedures developed and successfully applied for the (single impurity) 
Anderson model. We present a systematic study of the magnetic susceptibility 
and compare our results with those of a Hubbard model at large U . Similarities 
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and differences are pointed out and tfie magnetic phase diagram of the t-J model 



is discussed. 



Pacs numbers: 71.27+a, 71.28+d, 75.10.LP 
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1 Introduction and survey. 

The description of strongly correlated electron systems involves by and large three 
different classes of models. First one may consider a system consisting of uncorre- 
lated delocalized electronic states hybridizing with localized states subject to a strong 
Coulomb repulsion. This situation is modeled by the well known periodic Anderson 
model frequently used to describe the so-called heavy- fermion compounds [|2[. The 
second important situation occurs when the delocalized states themselves feel locally 
such a strong repulsion. In that case one is led to the single-band Hubbard model 0, 
originally set up to describe (ferro-) magnetism and metal-insulator transitions in 3d 
transition-metals compounds like V2O3 but recently also used for the high-Tc super- 
conductors. Another interesting kind of system is obtained if in addition to those local 
correlations a nonlocal magnetic exchange is included. This is the domain of the so- 
called t-J model ^ which is frequently taken as an alternative to the Hubbard model 
to describe the properties of the cuprate superconductors. It is this model we want to 
study more closely in this paper. Although the t-J model may be viewed as an effective 
Hamiltonian for the low-energy properties of the Hubbard model in the limit of large 
local Coulomb energy i.e. vanishing effective magnetic exchange, both models are 
expected to differ fundamentally for increasing exchange interaction. 
The Hamiltonian of the t-J model reads 



In equation (|T]), X}^\^, = \j, M){j, M'\ are the standard Hubbard operators |^ acting 



t* 



(1) 



Ht-j 




3 



on states with quantum number M G {0, la} on site j, i.e. double occupancy of a 
site is strictly forbidden, and S** denotes the spin operator on site i. The sums in the 
Hamiltonian (0) are on nearest neighbors only. The transfer and exchange integrals t 
and J have been rescaled with the coordination number Z of the system to guarantee 
a physical meaningful result for large spatial dimensions to be introduced later. Note 
that for J* = the model (|I]) is the Hubbard model in the limit U = oo. An additional 
density-density interaction frequently included in the model (|l]) has been dropped here 
for reasons of convenience. 

Although the model (|l]) looks rather simple, relatively little is known exactly about 
its properties. In contrast to the Hubbard model, it is not even exactly solvable in d = 1 
except for the two special points J* = (Hubbard model) and J* = 2t* (supersym- 
metric t-J model) Nevertheless, exact diagonalization studies showed that the t-J 
model for (i = 1 and T = is a Luttinger liquid for all J < J^^, while for J > J^^ 
one finds phase separation into an electron and hole rich region Interestingly, close 
to this boundary, the ground state of the t-J model is dominated by superconducting 



pair correlations |]10|, while for smaller J antiferromagnetic correlations are strongest. 

Obviously, this would make the t-J model an interesting candidate for explaining 
e.g. high-temperature superconductivity. Unfortunately, the results for c? = 1 suggest 
a much too large value of J* /t* ~ 3 . . . 4 for this scenario. The interesting question 
thus is how these features survive in d > 1 and especially to what extent phase sepa- 
ration might occur at much lower values of J, as suggested by e.g. high-temperature 



expansions [ITTl . 



While in d = 1 the combination of exact diagonahzation and tools of conformal field 
theory provides a powerful framework to extract informations about the asymptotics of 
the macroscopic system, similar methods do not exist in d > 1. Quantum Monte Carlo 
techniques, too, cannot be applied for realistic lattice-sizes and temperatures due to 
a severe minus-sign problem. Thus most informations about the properties of the t-J 
model come from high-temperature expansions, which are restricted to relatively large 
values of J* and T [ll,!^, and exact diagonahzation studies for small two-dimensional 



systems [0,0 . The finite system size in the latter method possibly prevents one from 



resolving dynamically generated low-energy features, which one may especially expect 



close to half filling |jTj,|T5|,|T6| . Moreover, to interpret results for dynamic quantities 
calculated with this method one generally needs additional information from other 
techniques about the general structures to be expected. Clearly, a different approach 
to obtain results in the thermodynamic limit is needed. 

Usually, a mean-field theory provides a reliable tool to study at least the quali- 
tative features of models in theoretical solid-state physics. However, until recently a 
thermodynamically consistent mean-field theory like for spin systems did not exist for 
fermionic models like the t-J model (|l]): While the magnetic exchange term could in 
principle be handled by the standard Hartree factorization it is a priori not obvious how 
to treat the correlated hopping introduced by the first term in the model (HI) consis- 
tently within this ansatz. Different schemes, usually involving slave-boson techniques. 



have been proposed [Q. These methods treat the local dynamics induced by the corre- 
lations rather poorly and a systematic inclusion of fluctuations around the static limit 



to incorporate lifetime effects is very cumbersome and has not been successful yet . 

Over the past three years, however, a novel scheme was introduced to define a 
thermodynamically consistent mean-field theory for correlated systems that preserves 
the local dynamics exactly [^,^,^. In this contribution we shall use this so-called 
"dynamical mean-field theory" to study the mean-field magnetic properties of the t- 
J model (|l|). The paper is organized as follows. In the next section we will briefly 
introduce the dynamical mean-field theory and derive expressions for the magnetic 
susceptibility of the t-J model. We then present results on the magnetic properties and 
compare them to the large-?/ Hubbard model. A summary and discussion concludes 
the paper. 



2 Theoretical background 



Since the pioneering work of Metzner and VoUhardt [T^ and subsequently Miiller- 



Hartmann [0, Brandt and Mielsch and Janis |^ it is known that a correlated 
lattice model can be mapped onto an effective impurity system in the limit d ^ oo. 
This is one consequence of the important aspect of this limit, namely that the ir- 
reducible one-particle self energy is purely local [l^ll^ and a functional of the local 
propagator only [ pO| , pT|p^ , p^ ] . This property can be used to rewrite the lattice problem 
in such a way that one is left with the solution of an effective single-impurity Anderson 
model (SIAM), where the free bandstates are replaced by an effective medium obtained 
from the full problem with the site under consideration removed PD|,PT|,P^,|25| . The one- 
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particle Greens function or equivalently the one-particle self energy of the system are 
then given by the corresponding quantities of the effective single-site problem. We 
shall see later, that one can also calculate the two-particle correlation functions of the 
lattice system with the help of those of the effective SIAM. Note that this effective 
theory preserves the dynamics introduced by the local correlations and thus is still 
highly nontrivial since there does not exist a complete solution for the SIAM. However, 
there exist at least different numerical exact techniques like quantum Monte Carlo and 



controlled perturbational approximations to solve this local model [jT4|,^. All these 
methods can then in turn be used to provide a solution of correlated lattice models 
in the thermodynamical limit. This approach has become known as the dynamical 
mean-field theory. The name is based on the observations that (i) the limit d = oo pro- 
vides a canonical starting point for the construction of a thermodynamically consistent 
mean-field theory and (ii) in contrast to a standard mean-field theory (like e.g. the 
one for the Heisenberg model) one obtains a complex, frequency dependent function 
as molecular field due to the dynamical nature of the local Coulomb repulsion. Note 
that with the same arguments one also finds that the contribution to the one-particle 
self-energy due to interactions like the spin exchange in the model (|I]) is given by the 
corresponding Hartree diagram only and thus is also purely local and in addition static 



1^. The latter statement means that for d = oo the t-J model in the paramagnetic 
phase (i.e. when (5**) = 0) is identical to the Hubbard model with U = oo. Regarding 
the one-particle properties in this regime we thus expect the well known features of the 
Hubbard model |T^. The situation of course changes as soon as one has a transition 
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into a magnetic state which will be discussed elsewhere 



2.1 Susceptibility for the t-J model 

For our purposes it is convenient to represent the transverse spin susceptibility of the 
t-J model as 

P UJn ,Ulm 

where Xgi^^n, i^m] i^n) is the spatial Fourier transform of the particle-hole propagator 

/3 /3 /3 /3 



{TrCi^{n)4^{T3)Cj^{T2)cj^{n))tj . (3) 

In equations (H) and (|]) iuJn and ium denote Fermi Matsubara frequencies and iun a 
Bose Matusbara frequency. Quite generally, by introducing the irreducible two-particle 
self energy T]j{iu!n,iuJm','i^i), the particle-hole propagator can formally be written 
as 

(4) 

Here, xfj^i^^n] i^n) = —Gij{iuJn)Gji{iuJn + it^n) represents the unperturbed part of the 
particle- hole propagator and Gij{iuJn) the full one-particle Greens function of the sys- 
tem. 

Using standard techniques of field theory [^, one can express the irreducible 
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particle-hole self energy as functional derivative of the one-particle self energy with re- 
spect to the one-particle propagator. In combination with the observation, that within 
the dynamical mean-field theory (DMFT) (i) the one-particle self energy is purely local 
and (ii) the exchange term J* enters the one-particle self energy only on the Hartree 
level it follows that the two-particle self energy aquires the particularly simple form 

J* 

T\l{iun, iujp] ivi) = — —S\i^j\,n.N. + T'^^iuJn, iujp] m) . (5) 

The non-trivial second term is the irreducible particle-hole self energy for J* = 0, i.e. 
for the U = oo-Hubbard model. Note that within the DMFT this quantity is also 
purely local pO[] ! 



Inserting the result (^) into the expression (^) and transforming into g-space, we 
obtain as transverse magnetic susceptibility of the t-J model in the DMFT 

P p 

+ ^ xf\i^n; ii^n)T^\iu;n, iujp] il'n)Xq{i^p, i^m] il'n) ■ 
p p 

(6) 

J* 

In equation (R) denotes the Fourier transform of — — ^ij-ji n.N.- For the case of a 

simple hyper-cubic lattice one e.g. obtains J^= — -r X! '^'^^{li ' '^)- 

1=1 

The susceptibility (BI) contains as one contribution the susceptibility of the Hubbard 
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model in the limit U = oo given by [p2| 



P p 

(7) 

It is now straightforward to show that with the help of expression (|^) equation (j^) can 
be rewritten as 

XqiiuJn,iuJm;i'^n) = ^' {'i'^n, i^m] l^n) (8) 



Pi P p 

Performing the sums on n and m in equation finally leads to the appealing result 



Xfi^^n) = Xr'i^^n) 



1 - j^xn^< 



(9) 

'qX^ 

as expression for the magnetic susceptibility of the t-J model in the DMFT. Thus the 
major ingredient in the susceptibility of the t-J model is the corresponding quantity 
of the HM for U = oo. One should also note that the expression (P) is very similar to 
the standard RPA result 

Xqiii^n] f/ = 0) = Xq{^^n] t/ = 0, J = 0) [1 - Jg-x,-(^i/„; f/ = 0, J = 0)]-' (10) 

for the corresponding noninteracting system. Thus, as far as the DMFT for the t- 
J model is concerened, the susceptibility is formally obtained by simply replacing 
Xqiivn', ?7 = 0, J = 0) by Xq{i^n', U = oo, J = 0) in the RPA-formulas. Let us emphasize 
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that this correspondence holds only on a formal level: The physical situation described 
by (^) is of course fundamentally different from the one modeled by ([T0|)! 



2.2 The spin susceptibility of the Hubbard model 



As already mentioned, the dynamic spin susceptibility of the Hubbard model in real 
space is within the DMFT given by [p2| 



+ l3 J2 Xu\^^n] ih'n)^^K'i'^n, iuJp] iz/„) x/^^ (icUp, iuJm] ij^n) ■ 



Equation (|Tl|) obviously also holds for the local susceptibility, i.e. 



(12) 



with the same r^^(zu;„, itUp; iz/;) as in equation (11). Combining equations (0) and (12), 
the susceptibility can be expressed by the local susceptibility through a matrix equation 



X 



q,l 



eff 



-1 



i q,l 



1 pX ioc,i ■ r 



■V(o)-i _ y(o) 

Aq,l ^ loc,l t 



X 



loci 



(13) 



With the definition 



A; = Ai(iu;„) = \ Y.Xioc{i^n,i^m\wi) and the symmetry rela- 



tion xioc{i^n,i^^m',i^i) = Xioc{i^^m,i^n', —it^i) followiug from the definition (H) we can 
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formally perform the frequency sums in (^) to obtain 

Xr(^^0 = xr (^^0 + ^Af ■ • (14) 

P 1 ~ ;3 «oc,« ■ r 

as the final result for the magnetic susceptibility of the Hubbard model in the framework 
of the dynamical molecular field theory. 

It is important to note that until now no explicit reference to the value of U has 
been made, i.e. equation ( p^lf ) is valid for all U . The form ([I^) for the susceptibility 



of the HM is especially convenient for computational reasons, because the outer sums 
on Matsubara frequencies have been performed exactly. These can pose numerical 
problems because [Xq^i\nm decays at most like l/(nm) for large n, m and one has to 
care for the correct time ordering in the final sums (cf. equations and (pH])). Whereas 
for the inner sums the products occuring there lead to an asymptotic behaviour like at 
least ~ and thus a well defined sum. 

2.3 The local spin susceptibility 

The only unkown quantity in equation (0) is the local susceptibility Xioc{i^^n, i^m] i^-'n) 
defined by 

13 13 P 13 

/3 



xU^^n,iu:m-,ii^n) = ^ J J Jdr^ J dr, 6-"^-^^'-^'^ e-''^^^^'-^^^ 6-'"-^^'-^^'^ 



< T^Qt(r4)cJ(r3)Q|(r2)4(ri) > (15) 

Within the DMFT, this function is obtained from the corresponding quantity of an 
effective SIAM with the band electrons replaced by the effective medium of the DMFT. 
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For finite U, the most successful way to solve the effective single-site problem and 
calculate functions like (|15|) is by Quantum Monte Carlo techniques |22]. However, 



since we are interested in the limit [/ = (X) in the current context, this technique is not 
available. On the other hand, for U = oo the time-ordered perturbation theory p6| 
provides a natural and easy access to local quantities. In this method one expresses 
all local quantities through the resolvents Po(io-) (^) of the unoccupied (occupied) ionic 
states. Of course, this theory cannot be solved exactly, so further approximations have 
to be introduced. Here, we shall use the so-called non-crossing approximation (NCA) 
2^,p7[] to calculate these resolvents and express further local correlation functions of 



interest. In previous publications we have already shown that the NCA provides a 
reliable approximation scheme to calculate such local quantities |p!4| , p2| , p^ . Applying 



the standard diagrammatic rules of this perturbational technique in conjunction 
with the NCA we obtain 

1 f dz 
Zioc J, 27r? 

for the local susceptibility. In equation ([T6|) , Zioc = J2 § 4^^~^^Pm{z) denotes the local 

M C 

contribution to the partition function and the contour C surrounds all singularities of 
the integrands counterclockwise. 



3 Results 
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3.1 General remarks 



The expressions (|14D and ([16D in principle still allow for the calculation of the dynamical 
susceptibility. Unfortunately, the derivation of equation (0) utilizes the representation 
of all quantities in Matsubara-space, i.e. one would be left with the awkward task 
to analytically continue the results to the real axis. This nontrivial problem is left 
for a future publication ||28[. In this contribution we will concentrate on the static 



susceptibility, i.e. we set iui = 0. 

Before we turn to the actual results for the U = cxo-Hubbard and t-J model let 
us first briefly discuss the special limit (n) = 1. In this case the model (|I|) becomes 
equivalent to the Heisenberg model and it is a straightforward task to calculate the 
molecular field expression for the static susceptibility, which reads 

Comparing this expression with the result for the t-J model in equation d), one sees 
that obviously x^l^^ ^ f for (n) 1. On the other hand, (3/2 is also exactly the 
value we expect for the local susceptibility in this limit, i.e. Xiic^ ~^ f for {f^) ~^ 
1. From this it at once follows that the second part in equation (JT^) will become 
negligible for (n) close to half filling. On the one hand this offers a rather sensible test 
for the numerics involved in calculating the susceptibility for the HM. In addition it 
provides an interesting approximate ansatz for the susceptibility of the t-J model by 
setting Xif^^ ~ Xk,c^ in this limit. Note that this also allows for a simple approximate 
calculation of dynamics since Xk>c\^) much easier to obtain than x^^^(co') given by 
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(p^. The latter observation is especially interesting in the light of recent studies by 
Scalapino et al. who analyzed the dynamical susceptibility for the twodimensional t-J 
model obtained from exact diagonalization and found that it was rather well described 
by a form like @ with x^^^(co') replaced by some local quantity . 



3.2 The Hubbard model 

Let us start by discussing the Lindhardt function 

1 



xf = Y.Gu,i^u;^)G^,{^u^n) • (18) 

LU„,k 

While the whole derivation was completely independent of the actual lattice structure, 
we now have to specify the meaning of the fc-sum. We here choose a simple cubic lattice 
in d dimensions, i.e. the coordination number is Z = 2d, and take the limit d ^ oo 



to use the simplifications arising in this limit With t* = 1 as the unit of energy, 
one then obtains for the single-particle DOS the well-known Gaussian form po{e) = 
exp(— e^)/y^ ||19| and one can also evaluate the /c-sum in equation (|1^) analytically 



r|,|g] to yield 

oo 



Xq 



^ ujn {i^n + ^ - S(iu;„) - e) • {itOn + - S(iu;„) -e-r]^-e' ■ ^1 - if^ 



(19) 



In relation ([T9|) , rjg' = J2 cos{qi ■a)/d and Ti{z) is the one-particle self energy of the HM 

(=1 

for a given U > 0. Note that the external wave- vector q only enters via the function 
rj^ which basically describes surfaces of constant energy in the simple cubic Brillouin 
zone. For presentational reasons, we shall choose the special vector q = g(l, 1, 1,1,.. .) 
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and use the number q with < g < tt as label rather than —1 < rj^ < 1. 
0.5 

0.375 

^ 0.25 

0.125 

'^O 0.2 0.4 0.6 0.8 1 " 

q 

Figure 1: Lindhardt function for t/ = ,4, 7 and U = oo 

The Lindhardt function for the HM for four different values of ?7 = 0, 4, 7 and 
[/ = oo at a filling (n) = 0.95 and for a low temperature T = 1/30 is shown in Fig. |I[ 
Note the different scales for U = (right scale) and U = 4 ,7 and U = oo (left scale)! 
Without looking at the details it is thus clear, that the correlations induced by U 
strongly suppress this quantity. In addition one can observe a dramatic change in the 
g-dependence with increasing U. While for U = one has a strong peak at g = tt 
due to the nesting property of the simple-cubic Fermi surface close to half filling this 
feature is strongly suppressed by the damping introduced by the correlations for U = 4 
,7 and U = oo. In addition there occurs a cross-over from the maximum in xf^ being 
at g = TT for small [/ to g = for ?7 = oo. Note also that in contrast to U = the total 
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g-dependence is rather weak in the other cases. 

From the previous observation one may deduce two things: First, since for ?7 = oo 
there is no net magnetic exchange between neighbouring sites, we expect from the 



flatness of that also x^^^ will be relatively flat as a function of q. In addition, the 
fact that is maximal at g = suggests that x^*^ for t/ = 00 will be enhanced at 



g = rather than at g ^ vr as expected and observed for U < 00 p2 |. 
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Figure 2: Susceptibility of the HM at U = 00 as function of q and filling for 
two different temperatures T = 1/5 and T = 1/40. 

This behaviour can indeed be seen in Fig. |^, where we have plotted fo^' two 
different temperatures as function of g and doping 6=1 — (n). The susceptibility was 
normalized to its value at 5 = 0, i.e. to xf^i^ = 0) = (3/2. Note that we always find 
^HM^^ > 0) < x^^K^ = 0)- From the form (P) for the susceptibility of the t-J model 
it then at once follows that also x^J > 0) < x^J {5 = 0) for all values of J* and g. 



This should be compared with results from high-temperature expansions for d = 2 [11 
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which suggest a pronounced maximum in the uniform susceptibihty around 6 = 15% 
produced by spin fluctuations not included in the current mean-field treatment. 

Another interesting feature in Fig. Q is that in all cases the variation with q is 
comparatively weak, becoming somewhat stronger for lower temperatures and with 
increasing doping 6. We also observe a slight maximum at g = that becomes more 
pronounced for lower temperatures but interestingly weakens with decreasing doping 
for T fixed. This observation is substantiated by a look at the doping dependence 
of x^^^ in Fig. ^ for the local (circles), ferromagnetic q = (squares) and antifer- 
romagnetic q = tt (diamonds) susceptibility for an inverse temperature j3 = 30. It 
is interesting to note that the antiferromagnetic susceptibility of the HM at f/ = oo 
is always very close to the local one, which can be understood by the fact that due 
to the mapping of the HM onto an equivalent impurity model the local susceptibility 
already contains most of the (nearest-neighbour) antiferromagnetic correlations. Since 
for f/ = oo there is no additional net magnetic exchange the nonlocal corrections only 
give a small renormalization. In contrast to this the renormalizations for the ferro- 
magnetic susceptibility are comparatively strong and definitely tend to enhance this 
quantity above both the local and antiferromagnetic susceptibility. These results have 
to be interpreted in the light of Nagaoka's theorem where in the presence of one 
hole a ferromagnetic state for the background is favoured from a minimization of the 
hopping energy in the correlated system but not as a result of a direct magnetic cou- 
pling. Obviously, our results suggest that sizeable ferromagnetic correlations still exist 
for a finite number of holes. However, so far we do not find any hint towards a ferro- 
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Figure 3: Susceptibility for U = oo and (3 = 30 as function of 6. 

magnetic instability at low temperatures close to half filling. This is consistent with 
the conjecture that for bipartite lattices - like the simple hyper-cubic lattice studied 



here - the critical hole density for the Nagaoka state should be 5c = ||31 |. 



3.3 Results for the t-J model 

Inserting the results for the susceptibility of the HM at U = oo into equation (j^) 
we obtain the susceptibility for the t-J model as function of q and J* as shown in 
Fig. ^ for (n) = 0.95 and (3 = 30. The explicit exchange now obviously favours the 
antiferromagnetic point q = tt and eventually leads to an antiferromagnetically ordered 
state for J* > J* ^ 0.085 for this particular parameter set. 

The temperature dependence of 1/xaf ^ specific value of J* = 0.067 and three 
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Figure 4: Susceptibility of the t-J model as function of q for various values of 
J at a doping 6 = 5% and (3 = 30. 

dopings 6 = 2%, 6 = 9% and 6 = 15% is collected in Fig. ^. The full curve marks 
for comparison the case 6 = 0, where one has exactly 1/xaf — 2 ■ (T — J*/2). As 
expected for a mean-field theory, close to the antiferromagnetic transition one finds a 
behaviour 1/xaf = ~ '^N)/Ceff in all cases with decreasing Neel temperature T/v 
and decreasing effective Curie constant Cg// for increasing S (see e.g. inset to Fig. |^). 
It is quite noteworthy that close to half filling (i.e. for 6 = 2%) this linearity extends 
up to rather high temperatures. However, with increasing doping one eventually finds 
appreciable deviations from this linearity for temperatures well above T/v. Both T/v 
and Ceff vary roughly linear up to 15% doping. We would also like to point out that 
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T 

Figure 5: Inverse susceptibility of the t-J model as function ofT for J* = 0.067 
and three dopings 6 = 2%, 6 = 9% and 6 = 15%. Close to the phase transition 
one observes XApiT) = (T — Tiy)/Ceff as expected for a mean-field theory. 
Note that for 6^0 the linear behaviour is observed up to T = It*. The full 
line represents half filling, where Xaf — 2 ■ (T — J*/2). The inset shows the 
dependence of the Neel temperature Tjy and effective Curie constant C^ff on 6. 

up to a doping of 5 = 15% we do not observe any tendency towards incommensurate 
order. 

Witli tlie metliod outlined above we are now able to calculate the phase diagram 
Tn{6,J*) for the t-J model. The results for dopings 6 < 15% and J* < 0.12 are 
shown in Fig. |^. One observes the expected increase in the Neel temperature T/v with 
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Figure 6: Phase diagram T]\f {6, J) for the t-J model. The dashed lines represent 
(linear) extrapolations of the phase boundaries to T = 0. The corresponding 
values of Jc{S) behave like Jc{S) ~ 5^, as shown in the inset. 

increasing J* and a - for larger S roughly linear - decrease as function of S. We may 
use this approximate linearity of T/v((5) to extrapolate the curves T/v((5) for a given J* 
to T = 0. This procedure allows us to obtain an extrapolation for the phase diagram 
J* (5) of the t-J model at T = 0. The result is shown in the inset to Fig. We find that 
J* (5) behaves rather accurately like J* ~ 6^. The phase diagram in Fig. || should be 
compared to the DMFT results for the Hubbard model in the strong coupling limit . 
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In reference |]32| the authors calculate Tn{S, U) up to t/ = 7t*, which would correspond 



to J* ~ 0.14 for the t-J model. They also observe an almost linear dependence of T/v 
on the doping 5 for large values of U . However, although the value of T/v for 5 ^ 
and the observed linearity agrees quite well with our results, the depression of Tjv as 
function of 5 for the Hubbard model ai U = 7t* is much faster than in our Fig. ^ 
In addition one encounters a transition into an incommensurate state for 6 ^ 12% in 
the Hubbard model. Currently it is not clear whether these deviations - especially 
the lack of an incommensurate magnetic order for large doping - between the results 
for the large-[/ Hubbard model and the t-J model are real or due to the additional 
approximations introduced by using the NCA to solve the effective impurity problem. 
One should keep in mind, though, that for finite U respectively J* the Hubbard model 
and the t-J model are expected to show different physical behaviour: The mapping of 
the Hubbard model to an effective model with magnetic exchange generates in addition 
to the exchange term included in the t-J model also more complicated couplings, like 
for instance a three-site term which is also of the order J* and may give rise to quite 
important corrections in physical quantities [p3 |. 



Finally we should like to use the observation that close to half filling the suscep- 
tibility for the HM is relatively flat with respect to q and obtain an approximation 
for the dynamical spin structure factor S{q,uj) = Qmxq{uj)/{1 — e~^^) by assuming 
X^^^ijj) ~ Xiic^i^) ill equation (^. This approximation avoids the cumbersome cal- 
culation of the g-dependent susceptibility for finite frequencies. As an example the 
result for J* = 0.035, T = 1/30 and 6 = 5% is shown in Fig. |^. As expected, the 
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Figure 7: Approximate result for S^J{u) for 6 = 5%, T = 1/30 and J = 0.035. 

maximum in SV(^u!) is found at g = tt and u; = and the intensity decays very fast 
with increasing energy for all q. Since this quantity or its value at g = tt and u = can 
be measured by neutron scattering or NMR relaxation ||3^, it is definitely necessary 



to study the dependence on doping, temperature and J* more systematically. This is 
left for a future publication. 
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4 Summary and outlook 

We presented a theory and results for the magnetic properties of the t-J model in the 
framework of the dynamical mean-field theory, which treats both the correlated hopping 
of the fermionic degrees of freedom and the nonlocal exchange coupling between the 



spin degrees of freedom on the same footing. As has been pointed out this approach 
ensures a thermodynamically consistent description of the properties of the system and 
especially does not introduce artificial phase transitions like e.g. in slave-boson mean 
field theories. 

One in our opinion particularly interesting result is that the dynamical susceptibility 
of the t-J model can be expressed in an RPA-like fashion by the susceptibility of the 
Hubbard model at U = oo (cf. equation (|^)). In addition the latter can be split into a 
local part plus a g-dependend renormalization which for low doping turned out to be 
relatively small and only moderately varying with q. We find that in the case J* = 
(i.e. U = oo) the absence of an explicit magnetic exchange leads to Xg=n ~ Xii^ 
and an interesting enhancement of the ferromagnetic correlations. This is in contrast 
to the HM at finite U, where the effective magnetic exchange J ~ t'^/U leads to a 
strongly enhanced susceptibilty at g = vr and a suppression at g = instead. However, 
for the situation considered here - simple hypercubic lattice with nearest-neighbour 
hopping only - we did not observe a tendency towards a magnetic instability at g = 
for finite doping, in accordance with results obtained by other groups. The occurence 
of an enhanced ferromagnetic susceptibility for Hubbard model in the limit U = oo 
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nevertheless motivates a more detailed investigation of the mean-field properties of the 
Hubbard model in this particular limit for different lattice structures and longer-range 
hopping. 

A finite magnetic exchange J* again strongly enhances the antiferromagnetic sus- 
ceptibility. When one further increases J* one eventually encounters a phase tran- 
sition into an antiferromagnetic phase at a critical value J*{T,S). From our results 
of Xaf^{T,6) we extracted the phase diagram T/v(5, J*)- We found that T/v increases 
monotonically as function of J* and - for fixed J* - decreases monotonically as func- 
tion of S. For larger doping 5 we observed that the curves Ti^{5) for different but 
fixed values of J* are almost linear. This linearity agrees at least qualitatively with 
DMFT results for the Hubbard model at finite U, where one finds a crossover from 
standard weak-coupling behaviour in T/v (5) for small U to an almost linear variation 
for U >7t*. However, in contrast to our results one observes a much faster depression 
of T/v as function of S and in addition a transition into an incommensurate phase for 
large 6. Especially the latter feature was not reproduced in our calculations. The lin- 
earity of T/v (5) finally allowed us to extrapolate our data to obtain an approximation 
for the magnetic phase boundary of the t-J model at T = 0. 

The relatively weak dependence of the susceptibility of the HM on q was used to set 
up an approximation for the dynamical susceptibility by assuming x^^{^) ~ Xii^i^): 
thus giving to some extent a microscopic justification of the results in reference [ p9|] . 
Since in addition Xiic^i^) can be calculated fairly easy from the effective single-site 
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problem we were able to present results for the dynamical spin structure factor SV(^uj). 
The general expected features, i.e. sharp maximum at g = tt and a; = 0, a shift of the 
maximum to finite u for q < n and a fast decay as a; > 0, are well reproduced. There 
are of course several questions left. First of all one should check the assumption of a 
nearly g^independent x^'^\^) carefully for several values of doping and temperature. 
Second a systematic study of S^^i^uj) as function of doping and temperature is clearly 
needed. Another important issue not yet addressed concerns phase separation in the t-J 
model, which among other problems requires e.g. the evaluation of the compressibility 
in the antiferromagnetic phase. Work along this line is in progress. 

Finally, one should stress again that the results presented here were calculated with 
a generalized mean-field theory or equivalently for the limit d = oo. This obviously 
means that their applicability to e.g. the t-J model in c? = 2 or (i = 3 is unclear. 
From high-temperature expansions or exact diagonalizations for d = 2 one knows for 
example that the static homogenous susceptibility shows a nonmonotonic behaviour as 
function of 6, which may be attributed to fluctuations induced by the spin-flip term in 
the model (p. Since the DMFT neglects this type of processes it is not too surprising 
that in our results we always observe a monotonic decrease instead. We thus do expect 
that the predictions of the DMFT will be modified not only quantitatively but most 
likely also qualitatively, especially for two-dimensional systems. 
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